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Abstract. We examine a class of braneworld models in which the expanding universe 
encounters a "quiescent" future singularity. At a quiescent singularity, the energy 
density and pressure of the cosmic fluid as well as the Hubble parameter remain finite 
while all derivatives of the Hubble parameter diverge (i.e., i/, 77, etc. — > oo). Since 
the Kretschmann invariant diverges {Rikim.R^'^^™ oo) at the singularity, one expects 
quantum effects to play an important role as the quiescent singularity is approached. 
We explore the effects of vacuum polarization due to massless conformally coupled fields 
near the singularity and show that these can either cause the universe to recoUapse 
or, else, lead to a softer singularity at which H, H, and H remain finite while H 
and higher derivatives of the Hubble parameter diverge. An important aspect of 
the quiescent singularity is that it is encountered in regions of low density, which 
has obvious implications for a universe consisting of a cosmic web of high and low 
density regions — superclusters and voids. In addition to vacuum polarization, the 
effects of quantum particle production of non-conformal fields are also likely to be 
important. A preliminary examination shows that intense particle production can 
lead to an accelerating universe whose Hubble parameter shows oscillations about a 
constant value. 



PACS numbers: 04.20.Dw, 04.50.+h, 04.62. +v, 98.80.-k 



Quantum effects, soft singularities and the fate of the universe 



2 



1. Introduction 

Braneworld models give rise to interesting new physical effects jU E]- The well- 
known Dvali-Gabadadze-Porrati (DGP) model, for instance, can lead to an accelerating 
universe without the presence of either a cosmological constant or some other form 
of dark energy [3 . Generalizations of the DGP model can result in a phantom-like 
acceleration of the universe at late times 01 This class of models can also lead to 
new cosmological behaviour at intermediate redshifts: the loitering [Hj and mimicry [2j 
scenarios provide examples of cosmologies which are close to LCDM at late times but 
can show significant departure from LCDM-like expansion at 2; >, few. In both cases, 
the age of the high-redshift universe turns out to be larger than in LCDM, while the 
redshift of reionization is lower. Whether the universe has properties which are easier 
to explain within the braneworld context is an intriguing possibility demanding further 
exploration. In this paper, we examine yet another property of braneworld models 
which does not have parallels in general relativity: the possibility that the universe may 
encounter a quiescent future singularity as it expands [S]. 

It is well known that the rate of expansion of the universe and its ultimate fate 
depend on the system of equations governing evolution (general relativity, scalar-tensor 
theory, braneworld theory etc.) as well as on the form of matter and, particularly, on 
its equation of state. In general relativity (GR), if one assumes that matter satisfies the 
strong energy condition (SEC) p + 3p > 0, then, within a Friedmann-Robertson- Walker 
(FRW) setting, the evolution of the universe is strongly dependent upon the spatial 
curvature: a spatially closed universe turns around and collapses whereas open and 
flat cosmologies continue to expand forever. The situation becomes more complicated 
(and interesting!) if one of the following assumptions is made: (i) the expansion of the 
universe is not governed by GR, (ii) matter can violate the SEC and even the weak energy 
condition (WEC) p + p > 0. In the latter case, if w = p/p < —1, then the expanding 
universe can encounter a "big rip" future singularity at which the density, pressure, 
and Hubble parameter diverge. In the former case, if the equations of motion have 
been derived from a braneworld action jl], then the expanding universe can encounter 
a different kind of (quiescent) future singularity, at which the density, pressure and 
Hubble parameter remain finite, but derivatives of the Hubble parameter, including 
H, diverge as the singularity is approached jH]. The occurrence of this singularity is 
related to the fact that the equations of motion are no longer quasi-linear (as they are, 
for instance, in GR) but include terms which are non-linear in the highest derivative. 
Examples of such singularities can be found in models other than the braneworld model; 
for instance, [Oj refer to singularities in which the deceleration parameter tends to infinity 
as the "Big Brake" while, in [TU], they are called "sudden" singularities (see also [TT]). 
The geometrical reason of a quiescent singularity in braneworld models is connected 
with the fact that the brane embedding in the bulk becomes singular at some point 
(see |B] for details). Since quiescent singularities can occur both in the past and in 
the future, they might provide an interesting alternative to the more conventional "big 
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bang" / "big crunch" singularities of general relativity. An important distinction between 
the quiescent and the sudden singularity is the following: The existence of the quiescent 
singularity does not require matter with unusual properties, hence, both density and 
pressure remain finite near this singularity. For the sudden singularity, on the other 
hand, the pressure diverges as the singularity is approached, implying the presence of 
matter with exotic properties. 

In this paper, we examine the issue of how quantum effects might influence a 
braneworld which encounters a quiescent singularity during expansion. It is well known 
that quantum effects come into play when the space-time curvature becomes large, as 
happens, for instance, in the vicinity of a black hole or near the Big-Bang and Big- 
Crunch singularities of general relativity [131 • Since RikimR^^^^ ^ oo as one approaches 
a quiescent (sudden) singularity, one might expect quantum effects to become important 
in this case too (see, for instance, [Stll4i|). As we demonstrate in this paper, quantum 
corrections to the equations of motion at the semi-classical level result in several 
important changes in the evolution of the universe: due to the (local) effects of vacuum 
polarization, (i) the quiescent singularity changes its form and becomes a much weaker 
"soft" singularity, at which H and H remain finite but i7 — > oo; (ii) vacuum polarization 
effects can also cause a spatially fiat universe to turn around and collapse. Both (i) and 
(ii) demonstrate that the incorporation of quantum effects into the braneworld equations 
of motion can radically alter the future of the universe and lead to behaviour which 
differs significantly from general-relativistic cosmology. Since the quiescent singularity 
(and the associated quantum effects) arise as the density drops below a threshold value, 
it follows that those regions which are significantly underdense (voids) may be the first to 
encounter the quiescent singularity. We also briefly discuss whether particle production 
effects could be significant as the universe approaches a quiescent future singularity. 
An interesting possibility that may arise in this case is that the universe expands in a 
regime in which the Hubble parameter vacillates about the de Sitter value. || 



2. Equations of motion 

We consider the simplest generic braneworld model with action of the form 

5 = [ / (7^ - 2Ab) - 2 / ^1 + / (m'R - 2a) + / L {hab, 0) . (1) 

.v/bulk v/brane J v/brane ^ ' v/brane 

Here, TZ is the scalar curvature of the metric Qab in the five-dimensional bulk, and R 
is the scalar curvature of the induced metric hat = gab — naUb on the brane, where 
is the vector field of the inner unit normal to the brane, which is assumed to 
be a boundary of the bulk space, and the notation and conventions of ^S] are used. 
The quantity K = h°'^Kab is the trace of the symmetric tensor of extrinsic curvature 

II It should be noted that braneworld models which approach the first quiescent singularity in the future 
(and have Friedmann-like behaviour in the past) appear to be in some tension with recent observational 
data However, the observational status of the oscillatory model discussed in Sec. [3 remains to be 
studied. 
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Kab = h^a^cnb of the brane. The symbol L{hab, 0) denotes the Lagrangian density of 
the four-dimensional matter fields (p whose dynamics is restricted to the brane so that 
they interact only with the induced metric hab- All integrations over the bulk and brane 
are taken with the corresponding natural volume elements. The symbols M and m 
denote the five- dimensional and four- dimensional Planck masses, respectively, Ab is the 
bulk cosmological constant, and a is the brane tension. 

Action leads to the Einstein equation with cosmological constant in the bulk: 

Gab + Kgah = , (2) 

while the field equation on the brane is 

m^Gab + (yhab = Tab + {K^b - habK) , (3) 

where Tab is the stress-energy tensor of matter on the brane stemming from the last 
term in action (jT)), i.e., 

1 5 / L{hab,(p) 

v/ brane 



- ab 



Vh Sh"-^ 

The five-dimensional bulk, satisfying Eq. is described by the metric 



(4) 



dsU = -f{r)dr^ + j^^+r'dQl, /W = - ^ - ^ . (5) 



Here, k = 0, ±1 is the sign of the spatial curvature of the brane, dflf. denotes the metric 
of the maximally symmetric three-dimensional Euclidean space with constant curvature 
corresponding to k, and the constant C, if it is nonzero, corresponds to the presence of 
a black hole in the bulk. The trajectory of the brane in the bulk is given by r = a(r), 
and then the brane is made the boundary by discarding either the region r > a(r) or 
the region r < a(r) from the bulk, resulting in two possible cosmological branches (see 
Eq. (Uni) below). 

The cosmological evolution on the brane that follows from Eqs. ((21) and © can be 
encoded in a single equation [H [TE] 



2 , K p + ay _ 4 , K Ab C 



where H = a/ a is the Hubble parameter, and p is the matter energy density on the 
brane. Here and below, the overdot derivative is taken with respect to the cosmological 
time t on the brane, which is connected with the bulk time r through the obvious 
relation 



dt 



{da/drf 
dr-\f^''^ f{a) ■ 

The term containing the constant C describes the so-called "dark radiation." The length 
scale i is defined as 

2m^ 
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Figure 1. Plot of relation in the case Ab > 0. Case (a) corresponds to pc > 0, 
and case (b) corresponds to pc < 0, where pc is given by Eq. H12|l . 



In what follows, we consider a spatially flat universe (k = 0) without dark radiation 
(C = 0). Then Eq. (jH)) takes the form 



3m^ 



4 ^/72 _ Ab 
6 



(9) 



This equation can be solved with respect to the total energy density on the brane 
Ptot = p + cr: 

2 



p + a 



1^2 



Ay 



(10) 



Sm^ £ V 6 

The "±" signs in the solution correspond to two branches defined by the two possible 
ways of bounding the Schwarzschild-(anti)-de Sitter bulk space by the brane, as 
described above [T3IIH1- Discarding the region r > a(r) (the region r < a(r)) from the 
bulk corresponds to the "+" ("— ") sign in (fTUj) . 

Alternatively, Eq. ^ can be solved with respect to H'^ with the result [H [13 IT9] 



1±\1 + 



p + cr _ Ab 
Sm^ 6 



Models with the lower ("— ") sign in this equation were called Brane 1, and models 
with the upper ( "+" ) sign were called Brane 2 in , and we refer to them in this way 
throughout this paper. 



3. Classical dynamics of the Brane 

Before we study quantum corrections to brane equations of motion, we describe possible 
classical dynamical regimes. Classical dynamics depends significantly on the sign of the 
bulk cosmological constant Ab. In the ensuing discussion, we shall examine separately 
all three cases, namely Ab > 0, Ab = 0, and Ab < 0. 
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• The case Ab > is shown in Fig.[T] The graph of © in the (if^, ptot) plane in Fig.[T] 
illustrates that in an expanding universe the matter density p decreases (except for 
a "phantom matter" which we do not consider in the present paper) , and the point 
in the plane (if^,ptot) moves from right to left in Fig. ^ 

A striking feature of Fig. ^ is that the value of the Hubble parameter in the 
braneworld can never drop to zero. In other words, the Friedmann asymptote 
if ^ is absent in our case. The upper and lower branches in Fig. ^ describe the 
two complementary braneworld models: branches AB and DB are associated with 
Brane2 and Brane 1 of |3], respectively, while branches AC and DC correspond 
to the lower and upper signs in (fTIH) . respectively, and describe the two branches 
with different embedding in the bulk. It should be noted that, in many important 
cases, the behaviour of the braneworld does not have any parallel in conventional 
Friedmannian dynamics (by this we mean standard GR in a FRW universe). For 
instance, the BC part of the evolutionary track corresponds to "phantom-like" 
cosmology with H > 0, even though matter on the brane never violates the weak 
energy condition. 

We would like to draw the reader's attention to the fact that, for a given ptot, a 
solution ii'^(ptot) exists if and only if ^ 



(see Fig. IH). This leads to two distinct possibilities for the late-time cosmological 
evolution of the braneworld. (i) If a > pc, then nothing prevents the matter density 
from diluting to p ^ at late times. In this case, the braneworld approaches a 
De Sitter- like future attractor at which H —>■ constant, (ii) In the opposite case, 
when a < pc, the braneworld dynamics is very different. In this case, Eq. (fT^ 
can be rewritten as p > pc — o", which implies that the density of matter can only 
drop to Pc — cr and no further ! Indeed, at p —* — cr, the universe experiences a 
"quiescent" singularity, at which the density p and the Hubble parameter H remain 
finite, while derivatives of the H, including H, H etc., diverge [Hj. (For Ab < 6£~^, 
we have pc < 0, and, in order for the quiescent future singularity to exist, the brane 
tension a must be negative (see Fig. [TJd).) 

• The case Ab = is shown in Fig |21 The Minkowski bulk (Ab = 0) leads to two 
generic future attractors, namely, the De Sitter-like Brane 2 jH] and the Friedmann- 
like braneworld evolving to the point C (with H —>■ and p 0). However, 
a sufficiently large negative brane tension (cr < — 3m^/£^) can also trigger the 
formation of a quiescent future singularity at the point B. 

• The case Ab < is shown in Fig. El A very interesting situation arises in this 
case since the Hubble parameter can reach zero within a finite interval of time with 
fl < subsequently. This situation describes the recoUapse of the universe. 

It is worth stressing that the recollapse of the braneworld is a consequence of the 
modified expansion law and is not due to the presence of a spatial curvature 

% This corresponds to the quantity under the square root in Hll|l being non-negative. 




(12) 
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Figure 2. Plot of relation © in the case Ab = 0. The point C dividing the branches 
with different embedding corresponding to different signs in I|1U|) is at the origin of 
coordinates. 




Figure 3. Plot of relation Q in the case Ab < 0. Case (a) corresponds to the 
condition |Ab| < 6£^^, and case (b) corresponds to the condition |Ab| > 6^^^. The 
dashed part corresponds to the forbidden nonphysical branch. 

term — the recollapsing braneworld is spatially flat ! 

It is easy to show that Brane 1 (lower branch) invariably recoUapses provided the 
brane tension is less then po = M^-y/— 3Ab/2 , which is the energy density at the 
point D in Fig. The critical case o = pd leads to a Friedmann late-time 
dynamics [201 • The Brane 2 case (the upper branch AB), however, leads to two 
different possibilities. Furthermore, this difference is important if we allow for the 
existence of negative brane tensions. For instance, in Fig. a Brane 2 evolves 
along the AB part and encounters a quiescent singularity at the point B. This 
requires negative tension a < 3m^(Ab/6 — £~^) [H]. The part BC corresponds to 
a hypothetical universe which is born at the singularity B, expands to the point 
C, then contracts and ends up in the quiescent singularity at B. The BC track 
is characterized by the negative derivative of ptot as a function of H^. Simple 
calculations lead to the condition |Ab| < Qi~^ for the existence of this branch. For 
|Ab| > 6^~^, the situation is presented in Fig. Eb- In this case, the Brane 2 universe 
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also recoUapses (provided the brane tension is negative and sufficiently large by 
absolute value). 

4. Quantum corrections to the equations of motion 

It is well known that quantum effects, including vacuum polarization and particle 
production, generically occur within regions of strong space-time curvature such as in 
the vicinity of black holes and near the cosmological Big-Bang singularity. As we have 
seen earlier, the braneworld — in addition to having the usual Big-Bang and Big-Crunch 
singularities — also possesses quiescent singularities, which an observer can encounter 
while the universe is still expanding. This singularity is specific to the braneworld, 
since the density of matter, its pressure, and the Hubble parameter all freeze to 
constant values, whereas if, H etc. diverge as the singularity is approached. Since 
this particular singularity can also develop in an accelerating universe akin to ours [H], 
it is of considerable interest to ask whether the nature of the quiescent singularity will 
in any way be affected once quantum effects are incorporated into the treatment. 

In general, quantum effects in curved space-time can arise on account of the vacuum 
polarization as well as particle production. It is well known that the latter is absent 
for conformally invariant fields (which we shall consider in this section) and that, in 
this case, quantum corrections to the equations of motion are fully described by the 
renormalized vacuum energy-momentum tensor which has the form jTHj 

(T,,) = k.HlS + k.H'f^ + k.H''^ , (13) 

where 

H^^a = 2DaDtR - 2habD'D,R + h^bR - 2RRab , (14) 

2 1 1 

h!^^ ^ - R\R^b + -RRab + -habR^'Rc, - -habR' , (15) 

and if^^"* is a local non-geometric tensor which depends upon the choice of the vacuum 
state and has vanishing trace (ii^4)" = O); ki,k2,ks depend upon the spin weights of 
the different fields contributing to the vacuum polarization. 

Equations ()14|) and ()15|) lead to the following vacuum expectation value for the 
energy density: 

= (Too) = k2H^ + ks (2HH + QHH^ - H^) . (16) 

In order to assess the effects of the vacuum polarization on the dynamics of the 
braneworld, one must add pq to the matter density in (fTT|) or (fTT1|) so that p p + pq 
in those equations. An important consequence of this operation is that the form of the 
equation of motion changes dramatically — the original algebraic equation changes to a 
differential equation ! The dynamical equation (fTUI) now takes the form 

HH = ^H^ - 3HH^ + {2k^)-^ (^-k2H^ + 3m^H^ - ptot ± SM^if^ _ Ab/e) . (17) 
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The two signs in (jl7p correspond to two different dynamical equations. As explained 
in Sec. |21 the sign is fixed by specifying one of the two possible ways of embedding the 
brane in the bulk. 

The goal of the present paper is to study the stability of the classical solutions 
when vacuum polarization terms are taken into account. The /c2-term in which 
does not contain time derivatives of H, can only change the position of the future stable 
points. On the contrary, due to the k^-ieim. in fll6|) . some classical solutions can lose 
stability. Therefore, for simplicity (and without loss of generality), we set = in our 
calculations. This assumption simplifies the situation significantly since, in this case, 
the stationary points of Eq. (fTTjl are precisely the solutions of the classical equation © 
with p = 0. 

For m = (the Randall-Sundrum brane), Eq. |T7j) becomes 
1 



HH = -H^ - + [-k^H^ - ptot ± 3MY^' - Ab/6j / 2k,, (18) 

We can see immediately that the stationary points can exist only for the "+" sign in 
fll8|) . This statement is valid also for > 0. A negative value of ^2 leads to a De Sitter 
solution, which has no analogs in the classical case, so we do not consider it here. By 
linearizing the brane equation of motion near the classical stationary point, it is easy 
to find that the condition for this stationary point to be stable is ^3 < 0, which is the 
same condition as in the standard cosmology. 

We now turn to the case m 7^ 0. If the brane has nonzero tension a, the stationary 
points of (|T7jl in the case /c2 = can be found by substituting a into (jHl) and setting p = 0. 
After that, we linearize Eq. ()17|) at these stationary points and find the eigenvalues of 
the corresponding linearized system. The condition of stability of the stationary point 
is that its eigenvalues are negative. 

The eigenvalues at the stationary points where H = are given by 

pi,2 = l{fi± + 4/2) , (19) 
where we have made the notation 

h = -SH, (20) 

A--|i + 77;±— . (21) 



Two different signs in Eq. !^21\i correspond to two different equations of motion, 
while, in Eq. (|19|). we have two different eigenvalues of a single equation. 

Since /i is negative, the eigenvalue p2 corresponding to the "— " sign in Eq. (fT^ is 
also always negative. Moreover, pi is positive if and only if /2 is positive. As a result, 
the stability of a fixed point is equivalent to the condition /2 < 0. 

From now on, we consider only the second eigenvalue, p2, so the two signs in the 
expressions given below always corresponds to two different equations of motion. 



Quantum effects, soft singularities and the fate of the universe 



10 



The equations for the fixed points are 



(23) 



Substituting the value of A obtained from this equation into ()21|). we can rewrite it 
in the form 

2ri 



/2 




2 ± 



(24) 



'H^ - Ab/6^ 

We classify solutions corresponding to the upper sign ("+") and the lower sign 
( "— " ) as the "+" branch and "— " branch, respectively. Then, for the "+" branch, the 
expression in the brackets is positive, so we make our first conclusion: 

• Fixed points of the "+" branch are stable with respect to quantum corrections if 

^3 <0. 



Proceeding further, we solve Eq. ()23|) with respect to H . The solutions are given 

by 

= X + 2£-2 ± 2ryA - Ab/6 + £-2, (25) 

however, it is necessary to be careful about the two different branches. There are two 
possible cases: 

(i) A > Ab/6. In this case, we have one solution for the "+" branch (with the upper 
sign in (j25p ). and one solution for the "— " branch (with the lower sign in (j25|) ). 

(ii) Ab/6 — £^2 < A < Ab/6. In this case, two solutions belong to the "— " branch, 
while the "+" branch does not have fixed points. 

For A < Ab/6 — £~2^ there are no fixed points for both equations of motion. 

Consider the first situation. We already know that stability of a fixed point for the 
"+" branch requires < 0. For the "— " branch, we have 



/2 



1 

2h 



( 



2 - 



2^ 



-1 



V 



A 



(26) 



Ab/6 + 2£-2 + 2£-yA - Ab/6 + £-2 j 

at a fixed point. As we have now A — Ab/6 > 0, the second term in the brackets is always 
smaller than \/2, so the whole expression in the brackets is positive. This implies the 
following conclusion: 

• If the "— " branch has only one stable point, its stability requires < 0. 

For the second case, both solutions belong to the "— " branch. The upper sign 
gives the same expression for /2 as in ^I^. The condition A — Ab/6 + £~2 > g makes 
the second term in the brackets of (j^^ smaller then 2, and again the total number in 
the brackets is positive. On the other hand, for the lower sign in (j2SI), we have 



2t 



\ 



Ab/6 + 2£-2 - 2£-y A - Ab/6 + £-2 j 



(27) 
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One can see that the number in the brackets is negative. As a result, we arrive at the 
following conclusion: 

• If two stationary points belong to the "— " branch, then one of them is stable 
and the other one is unstable. 

Having in mind the description of the fixed points given in Sec. OJ we obtain the 
following picture: 

The simplest case corresponds to the classical points described by Fig. Eb- Note 
that now the horizontal axis represents a, because we are dealing with the stationary 
points only. Each equation, with "+" and "— " sign, has only one possible stationary 
point, and these are stable for /C3 < 0. 

For the range of the parameter Ab corresponding to Fig. EK, stationary points of 
the "— " sign of Eq. (fTTj) are located in the DE part. Their stability also requires k^, < 0. 
However, stationary points of the equations of motion with "+" sign are located in the 
whole admissible region AC. For a < pc (see Fig. 3a), we have two possible stationary 
points, one on the AB branch, and the other on the BC branch. Their stability properties 
are opposite to each other. In particular, the unusual condition ks > is required for the 
points on BC branch to be stable. This conclusion is, however, not so important because 
the branch BC cannot be reached as a result of conventional cosmological evolution in 
the classical picture. 

More interesting situation is realized for Ab > (see Fig. [T]). Equation (fTTj) with 
"— " sign again has stationary points on the CD branch, and that with "+" sign has 
it on the AC branch. If k^ < 0, points on the AB branch are stable, and those on the 
BC branch are unstable; for k^ > 0, the situation is opposite. However, now classical 
and quantum dynamics differs significantly. In the classical picture, the BC branch is 
part of the evolution track of the Brane 2, and it can be reached during a cosmological 
evolution from some high-energy initial phase. The effects of vacuum polarization make 
this branch unreachable for a braneworld with k^ < 0. 

These results can be interpreted from a different point of view. It is reasonable to 
assume that the evolution of a brane should be close to the classical picture in regions 
which are far removed from singularities. That is why we are interested only in the case 
/cs < 0. What is the fate of a braneworld if we take quantum corrections into account? In 
general, this problem requires numerical integration of Eq. ()17j) : however, we can make 
several qualitative statements. Fig. contains only stable branches, so the dynamics 
with quantum corrections is qualitatively the same, with recoUapse somewhere in the 
neighbourhood of the points B and D. In the case shown in Fig.lHt, the Brane 1 evolution 
DE does not change much and results in a recoUapse. The Brane 2 branch AB in the 
classical picture meets a quiescent future singularity with H finite and H —00 at 
the point B. Such a singularity is absent in the quantum picture because the expression 
under the square root in ()17|1 is positive if Ab < 0. Numerical simulations show that the 
ultimate fate of the Brane 2 universe in the case with Ab < and with a large negative 
a is also a recoUapse. 
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On the other hand, in the case Ab > 0, recollapse becomes impossible since the 
existence of a square root in ()17p requires > Ab/6. This restriction is vahd also 
for the classical equation (jH)). However, during the classical evolution, a braneworld 
reaches the point C, where = Ab/6, with H = 0, and then enters the CB branch, 
which corresponds to super-exponential expansion {H > 0) but which is impossible in 
the quantum case. Our numerical results show that, instead, the universe reaches some 
point with H'^ = Ab/6 and H <0. This causes the square root in ()17p to vanish, which, 
in turn, leads to the divergence of H, while H remains finite. ^ As a result, instead of 
a quiescent future singularity with \H\ oo at the point B, the universe encounters a 
much weaker singularity with ^ oo in the neighborhood of the point C. Similarly, 
the Brane2 universe misses the point B and falls into a singularity with H —>■ oo. 



5. Quiescent singularities in an inhomogeneous universe: a preliminary 
analysis 

The preceding discussion focussed on a homogeneous and isotropic universe whose 
expansion was governed by the brane equations of motion. Since the real universe 
is quite inhomogeneous on spatial scales < 100 Mpc, it is worthwhile to ask whether 
any of our previous results may be generalized to this case. 

Although we are not yet able to provide a self-consistent treatment of the brane 
equations for this important case, still, some aspects of the problem can be discussed at 
the phenomenological level. Consider, for instance, the expansion law ()11|) 



P 



a ^2 
3m^ £^ 



1±\ 1+P 



p + cr 
3m^ 



Ab 
6 



Ab 
6 



1 




(2^ 



A necessary condition for the existence of a quiescent singularity is that the matter 
density p drop to a value which is small enough for the square root on the right- 
hand side of ()28|) to vanish. When this happens, the universe encounters the quiescent 
singularity at which p and H remain finite, but d and higher derivatives of the scale 
factor diverge. Note, however, that, according to ^Ei), the universe encounters the 
quiescent singularity homogeneously, i.e., every part of the (spatially infinite) universe 
encounters the singularity at one and the same instant of time. This follows from the 
fact that the density in depends only upon the cosmic time and upon nothing else. 
In practice, however, the universe is anything but homogeneous, its density varying from 
place to place. For instance, it is well known that the density of matter in galaxies is 
> 10^ times the average value while, in voids, it drops to only a small fraction of the 
average value. This immediately suggests that the brane should encounter the quiescent 

+ A similar result was obtained in jl4| for the case of the Randall-Sundrum braneworld model 
containing matter with unusual properties resulting in sudden singularity. 
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singularity in a very inhomogeneous fashion. Underdense regions (voids) will be the first 
to encounter the singularity. Even in this case, since the density in individual voids is 
inhomogeneously distributed, more underdense regions lying closer to the void center 
will be the first to experience the singularity. It therefore follows that the quiescent 
singularity will first form near the centers of very underdense regions. As the void 
expands, its density at larger radii will drop below p^, where 



consequently, the singularity will propagate outward from the void center in the form 
of a quasi-spherical singular front. (For simplicity, we have assumed that all voids have 
a spherical density profile; this assumption may need to be modified for more realistic 
cases; see j2Il 122] and references therein.) 

The above approach provides us with a very different perspective of the quiescent 
singularity than that adopted in the previous sections (and in [H]). For one thing, the 
singularity may be present in certain regions of the universe right now, so it concerns 
us directly (as astrophysicists) and not as some abstract point to which we may (or 
may not) evolve in the distant future. The second issue is related to the first, since the 
singularity could already exist within several voids (there are as many as a million voids 
in the visible universe in at least some of which the condition p — Ps could be satisfied), 
a practical observational strategy needs to be adopted to search for singularities in 
voids. (Similar strategies combined with strenuous observational efforts have led to the 
discovery of dozens of black holes in the centers of galaxies |23j.) 

A number of important issues therefore need to be addressed: 

(i) Since RikimR^^^^ ~^ oo within a finite region at the very center of a void, it follows 
that, unless this region is contained within an event horizon, we will find ourselves 
staring at a naked singularity ! (As shown earlier, quantum effects do soften the 
singularity so that RikimR^^''"^ may remain finite if these effects are included.) 

(ii) In the discussion in Sec. HI the issue of particle production was ignored since it 
was assumed that we were dealing with conformally coupled fields which are not 
created quantum mechanically in the (conformally flat) homogeneous and isotropic 
universe which we have been considering up to now. However, the moment we drop 
the homogeneity assumption, the issue of particle production immediately crops 
up, and we must take it into account if our treatment is to be at all complete 
|24j . (In a related context, the quantum creation of gravitons takes place even in 
a homogeneous and isotropic universe, since these fields couple minimally, and not 
conformally, to gravity |25j.) 

Let us discuss the possible effect of particle production in more detail. First, 
we consider the model of homogeneous universe taking it as an approximation to 
the situation inside an underdensity region (void). Gravitational quantum particle 
production occurs as the singularity is approached. Since the local value of the Hubble 




(29) 
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parameter remains finite at the singularity, production of tlie conformally coupled 
particles (like photons) is expected to be negligible. However, particles that are non- 
conformally coupled to gravity (which could be, for example, Higgs bosons in the 
Standard Model) will be copiously produced as the acceleration of the universe a rapidly 
increases. The rate of particle production depends not only on their coupling to gravity 
but also on their coupling between themselves. Gravitationally created primary particles 
will decay into conformally coupled secondaries (electrons, photons, neutrino, etc.), 
which will influence the rate of production of the primaries by causing decoherence 
in their quantum state. The whole process is thus not easy to calculate in detail. 
However, from very general arguments it can be seen that creation of matter due to 
quantum particle production is important for the dynamics of the universe during its 
later stages. * 

For the sake of physical simplicity, we restrict ourselves to the case of vanishing 
bulk cosmological constant Ab and write Eq. ()28|) in the form 




(30) 



3^2 

Ap = p- ps, ps = -(T- , (31) 



g= :r (32) 



where 



and where we have chosen the physically interesting "+" sign in Eq. (j28|) . Thus, we 
have two free parameters in our theory, namely, I and ps- The value of m is assumed to 
be of the order of the Planck mass. In this case the early-time behaviour of the universe 
follows the standard Friedmann model, as can be seen from (j28|) or (j3(J|) . 

Let the average particle energy density production rate be Pprod- Then, 
differentiating Eq. (jHUI) . we obtain 

Pprod - iH p 
2{?,m?/\pf 

where 7 > corresponds to the effective equation of state of matter in the universe: if 
p = wp, then 7 = 3(1 + w). The second term in the numerator of ()32|1 follows from 
the conservation law and describes the effect of the universe expansion on the matter 
density. (Note that p includes contributions from quantum and classical matter.) 

In order to qualitatively assess the effects of particle production, let us examine 
two fundamentally distinct possibilities. 

(i) Suppose that, in the course of evolution, Ap — is reached in a finite interval of 
time. Since the Hubble parameter is a unique function of the energy density, given 
by (jHUj) , and since the singularity value ps of the energy density is approached from 

* Effects of particle production are negligible in the neighbourhood of the usual cosmological singularity 
of the Friedmann universe because the energy density of ordinary matter and radiation strongly diverges 
and thus dominates at this sing ularity [211131 . In our case, the energy density of ordinary matter 
remains finite during the classical approach to the quiescent singularity, hence, particle production 
effects are of crucial significance. 
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above, it follows that i/ < in the neighbourhood of the singular point. In the 
purely classical case we find, after setting Pprod to zero in (jH^ . that H —oo as 
the quiescent singularity is approached. It is well known, however, that particle 
production effects are sensitive to the change in the rate of expansion jTSI, and it 
is expected that Pprod will go to infinity as if — > —oo. Since Pprod ^ jHp, this will 
result in H becoming positive, which contradicts the assumption that H < 0. 
Therefore, under the assumption that the critical density ps is reached in a finite 
time, the only possibility for H is to remain bounded. In other words, the rate of 
particle production should exactly balance the decrease in the matter density due 
to expansion, turning the numerator in ()H2|1 to zero: 

Pprod - iHp ^0 Pprod as p ^ Ps- (33) 

In this case, the universe reaches its singular state with the energy density due to 
particle production exactly balancing the density decrease caused by expansion, as 
given by (jH^ . 

(ii) It is not clear whether the above regime will be realized or whether, if realized, it 
will be stable, since it requires the exact balancing of rates (jH^ at the singularity. 
A second distinct possibility is that, due to the presence of particle production, the 
value of Ap = p ~ Ps always remains bounded from below by a nonzero density. 
In this scenario, H initially decreases {\H\ increases) under the influence of the 
increasing factor 1/ ^J Ap in (jH^ . However, a large value of \H\ induces active 
particle production from the vacuum which leads to an increase in the value of 
Pprod in ()32|). As the value of Ap reaches its (nonzero) minimum, we have if = at 
this point, according to (IHUj). after which the rate H becomes positive due to self- 
sustained particle production that continues because of the large value of the second 
time derivative H. After a period of extensive particle production, the universe 
reaches another turning point H = after which is continues to expand according 
to ()30p with decreasing energy density. Thus, we arrive at a model of cyclic 
evolution with periods of extensive particle production alternating with periods 
of classical expansion during which quantum particle production is negligible. This 
scenario bears a formal resemblance to quasi-steady-state cosmology proposed in 
a very different context by Hoyle, Burbidge, and Narlikar [2E]- The particle 
production rate in our case is estimated by the quantity Pprod given in (jHHj) . which 
is approximately the value it takes at the turning points where H = 0. The Hubble 
parameter in this scenario periodically varies being of the order of ii ~ £~^, and 
the energy density is of the order ps, so that particle production rate is 

Pprod ~ ^ • (34) 

Our discussion thus far was limited to quantum processes within a single underdense 
region (void) which was assumed for simplicity to be perfectly homogeneous. Let 
us now (qualitatively) discuss whether this scenario can be generalized to the real 
(inhomogeneous) universe. Clearly, the particle production rate Pprod in this case should 
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be regarded as being averaged over the entire universe, to which several significantly 
underdense voids are contributing. Equation (jH^ should therefore be treated as 
an ensemble average, where the mean particle production rate depends upon the 
distribution as well as dynamics of local underdensity regions. As a result, equation 
fjH^ is not expected to explicitly depend upon the behaviour of the Hubble parameter 
and, in principle, particle production can proceed even in a De Sitter-like universe, in 
which the Hubble parameter H remains constant in time. The rate of particle production 
in this case is given by equality with zero left-hand side: 

Pprod = iHp . (35) 

The value of the Hubble parameter in such a steady-state universe can be related 
to the f2-parameter in matter 

where we have used the basic Eq. (|30|l . For the average energy density, we obviously 
have p — ps ^ p. Hence, 




(37) 



or, finally. 



H ^ — — ^ . (3^ 



In principle, one might use these preliminary results to construct a braneworld 
version of steady-state cosmology, in which matter is being created at a steady rate in 
voids rather than in overdense regions (as hypothesized in the original version [2^]). This 
would then add one more model to the steadily growing list of dark-energy cosmologies 
12] • These conclusions must, however, be substantiated by a more detailed treatment 
which takes into account the joint effect of vacuum polarization and particle production 
near the quiescent singularity. Clearly, whether one or the other effect dominates will 
depend upon the number of conformal and non-conformal fields contributing to the 
vacuum, their spin weights, etc. We propose to return to this issue in a future work. 



6. Conclusions 



Cosmological models based on braneworld gravity have the interesting property of giving 
rise to singularities which are not commonly encountered in general relativity. This is 
largely due to the possibility of different kinds of embedding of the brane in the higher 
dimensional (bulk) space-time. Singular embedding implies that the expansion (in time) 
of the brane cannot be continued indefinitely. The singularities which we have examined 
in this paper arise because of this reason. They have the property that, while the density, 
pressure, and Hubble parameter on the brane remain finite, higher derivatives of the 

tt For comparison, the late-time value of the Hubble parameter in LCDM is 27 H = Hq\/1 — . 
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Hubble parameter blow up as the singularity is approached. For this reason, we refer 
to these singularities as being "quiescent." Despite its deceptively mild nature, the 
quiescent singularity is a real curvature singularity at which the Kretschmann invariant 
diverges {RikimR^^^^ oo). The importance of quantum effects in regions of large space- 
time curvature has been demonstrated in a number of papers jlH], and it should therefore 
come as no surprise that these effects can significantly alter the classical behaviour near 
the quiescent singularity, as demonstrated by us in this paper. Examining the vacuum 
polarization caused by massless conformally coupled fields, we find that these weaken 
the quiescent singularity resulting in an exceedingly mild (soft) singularity at which 
the values of iJ, H, and H remain finite while H and higher derivatives of the Hubble 
parameter diverge. 

Unlike the classical Big-Bang singularity, the quiescent singularity in braneworld 
models is reached in regions of low density and is therefore encountered during the 
course of the universe expansion rather than its collapse. Densities lower than the 
mean value are known to occupy a large filling fraction within the cosmic web I28j. 
Therefore, if the braneworld model is a true representation of reality, one might 
speculate that it is likely to encounter the quiescent singularity (or its quantum-corrected 
counterpart, the "soft singularity") within large underdense regions, or voids. The 
rapidly varying space-time geometry near the quiescent singularity can, in addition to 
vacuum polarization, also give rise to quantum creation of fields which do not couple 
conformally to gravity. (Examples include massive (leptons, higgs, etc.) as well as 
massless (gravitons) particles.) A preliminary estimate made by us in this paper shows 
that, if this process is sufficiently intensive, then the universe will expand at late times 
in a manner which is reminiscent of quasi-steady-state cosmology, with the Hubble 
parameter showing oscillations about a constant value. A more detailed estimate of 
particle production, however, lies beyond the scope of the present paper, and we hope 
to return to it in a future work. 
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